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OBJECT  OF  TASK 

To  Improve  exIsHng  knowledge  of  gamma  and  neutron  shielding  properties  of 
shelters  in  order  to  fill  gaps  in  nuclear  shielding  knowledge. 


ABSTRAQ 

This  study  investigates  the  additional  nuclear  shielding  from  an  isotropic  (plane) 
radioactive  gamma  source  afforded  by  various  shapes  of  underground  curved-roof 
shelters  compared  to  the  basic  slab  shield.  This  additional  shielding  is  defined  in 
the  form  of  a  dimensionless  Geometry  Factor  which  is  a  function  of  the  physical 
dimensions  and  shape  of  the  shelter. 

Curves  are  presented  from  which  the  Geometry  Factors  for  the  underground 
shelter  shapes  of  spheres^  horizontal  cylinders,  ellipsoids,  and  vertical  cylinders  or 
silos  can  be  obtained  with  minimum  calculations  using  only  the  physical  dimensions 
of  the  shelter  and  the  depth  of  material  above  the  crown  of  the  shelter. 

The  Geometry  Factors  thus  obtained  are  independent  of  material  and  photon 
energy  except  in  the  case  of  the  silo  and  can  be  used  as  dimensionless  factors  once 
the  attenuation  for  a  slab  shield  has  been  calculated  for  a  particular  set  of  radio¬ 
logical  conditions  and  materials. 
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Figure  1.  Skefch  of  most  general  case  considered.- 


INTRODUCTION 


The  U.  S.  Naval  Civil  Engineering  Laboratory  is  charged  by  the  Bureau  of 
Yards  and  Docks  with  the  responsibility  of  developing  design  principles  for  the 
construction  of  atomic  warfare  shelters  for  use  by  the  Naval  Shore  Establishment. 

The  simplest  shielding  configuration  for  nuclear  radiation  is  a  plane  slab,-  and 
this  case  has  been  taken  as  basic.  Much  work  has  been  done  for  slabs  with  both 
isotropic  and  plane  collimated  sources  at  various  angles  of  incidence.  Since  actual 
underground  shelter  geometry  may  be  quite  different  from  slab  geometry/  the  shield¬ 
ing  effectiveness  of  other  configurations  needs  to  be  known.  This  report  is  a  study 
of  the  effectiveness  of  underground  structures  of  various  curved  roof  shapes  as  shields. 


BASIC  ASSUMPTIONS 

The  following  assumptions  have  been  made; 

1.  The  shelter  is  buried  in  a  semi-infinite  mass  of  earth  so  that  the  distance 
from  the  crown  of  the  structure  to  the  surface  of  the  ground  is  t^^. 

2.  The  source  of  radiation  is  a  plane,  equally  distributed,  isotropic  source  of 
gamma  photons  located  on  the  surface  of  the  ground. 

3.  The  surface  of  the  ground  is  smooth  and  infinite  in  extent. 

DESCRIPTION  OF  THE  PROBLEM 

It  is  assumed  that  there  is  a  shelter  buried  below  an  infinite  plane  source  of 
gamma  photons.  Figure  1  is  a  sketch  of  the  most  general  case  considered.  This 
shape  is  a  hemi-ellipsoid  whose  axis  radii  are  h,  b,  and  c.  In  all  cases  only  the 
dose  to  the  central  point,  P,  in  the  shelter  will  be  considered.  One  method  to 
determine  the  total  dose  to  point  P  from  the  entire  infinite  plane  is  to  sum  up  the 
contribution  from  all  points  in  the  plane. 
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Referring  fo  Figure  2,  which  is  fhe  plane  section  through  P  and  P^/  it  can  be 
seen  that  the  shield  mass  con  be  considered  in  two  parts:  (1)  a  slab  shield  of  thick¬ 
ness  t^  and,  (2)  an  additional  mass  bounded  by  the  slab  shield  and  the  surface  of  the 
shelter.  The  distance  t  of  the  exponential  term  of  Equation  2  can  be  separated  into 
two  terms,  t-j  and  f2»  Regardless  of  the  shape  of  the  shelter,  as  long  as  it  is  concave 
downward,  the  value  of  t-j  can  be  expressed  as  a  function  of  t  and  the  central  angle 
0  since  it  is  symmetrical  about  the  vertical  axis,  thus: 


(3) 


By  similar  reasoning,  the  distance  from  the  dose  point,  P,  to  the  underside  of 
the  slab  shield  would  be:  h  sec0.  Converting  to  spherical  coordinates,  the  distance 
R^  (dose  point  to  roof  of  shelter)  would  be: 


R 


1 


2  2 
^  sin  0  cos  Q 


9  9 

sin^  0  sin^  Q 


cos^  0 


-1/2 


(4) 


the  distance  t  in  Equation  2  can  now  be  expressed  as: 


t  =  t  sec  0  +  h  sec  0  -  Ri 
m  1 


(5) 
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Figure  2.  Plane  secHon  through  points  P  and  ?2  of  Figure  1. 
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Since  only  concave  downward  cases  will  be  considered  in  fhis  report',  the 
special  case  of  the  uniform  slab  shield  of  infinite  extent  will  provide  the  least 
attenuation  of  the  gamma  radiation.  The  slab  shield  then  will  be  the  basic  configu¬ 
ration  against  which  all  other  cases  will  be  compared.  The  geometry  factor,  GF, 
will  be  defined  as: 


G  F  = 


where:  Dgj  =  gamma  dose  through  shield  of  geometric  shape  being  considered 

Dj  =  gamma  dose  through  slob  shield  of  thickness  t^ 

Thus,  only  the  factor  under  the  integral  sign  of  Equation  2  need  be  considered 
further,  since  Sf/2  will  be  common  to  all  cases. 


BUILD-UP  FAaOR 


The  build-up  factor  is  the  most  questionable  quantity  in  the  numerical  computa¬ 
tion  of  the  geometry  shielding  factors  for  shelters  of  various  shapes.  It  has  been 
defined  as  the  ratio  of  some  measurable  property  of  the  photon  beam  (i.e.,  inten¬ 
sity,  number  of  photons,  energy-flux,  or  biological  dose),  when  the  effects  of  all 
quanta  are  included,  to  the  measurable  property  obtained  when  only  the  uncollided 
flux  is  considered.  Consequently,  it  is  possible  for  the  sake  of  mathematical  devel¬ 
opment  to  define  the  build-up  factor  in  terms  of  the  geometry  being  considered.  A 
mathematical  formula  has  practical  application  only  when  numbers  are  substituted 
in  it  to  determine  some  real  quantity.  It  was  necessary  therefore  to  utilize  available 
build-up  factors  for  the  numerical  integration  used  In  this  report. 

To  arrive  at  the  attenuation  integral,  an  integration  or  numerical  summation 
over  point  Isotropic  sources  was  used.  Consequently,  the  dose  build-up  factors  for 
point  isotropic  sources  derived  by  the  moments  method  were  used.  ^  These  build-up 
factors  were  derived  from  the  assumption  that  the  source  and  detector  were  located 
in  an  infinite,  homogeneous  medium,  and  the  error  involved  In  applying  these 
build-up  factors  to  finite  and  even  specific  geometries  must  be  considered. 


in  fhe  actual  situation  being  considered,  where  the  detector  is  located  deep 
in  a  shelter  cavity  at  some  depth  within  a  semi-infinite  medium,  using  point  iso¬ 
tropic  build-up  factors  for  an  infinite  medium  introduced  little  if  any  error.  This 
can  be  qualitatively  explained  as  follows.  The  semi-infinite  assumption  should  use 
a  smaller  build-up  factor  than  the  infinite  medium  since  we  are  neglecting  the 
photons  which  are  emitted  upward  from  the  source  and  then  are  scattered  back  down 
toward  the  detector.  But  since  these  photons  have  been  scattered  at  large  angles 
and  have  been  degraded  in  energy,  they  will  be  absorbed  sooner  than  those  photons 
emitted  downward  toward  the  detector,  and  their  contribution  to  the  dose  at  the 
detector  for  depths  of  two  mean  free  paths  (at  the  degraded  energy)  should  be  almost 
negligible.  The  greatest  error  would  result  v/hen  a  vacuum  existed  above  the  semi¬ 
infinite  medium.  In  an  air-soil  medium,  photons  will  also  be  back-scattered  from 
the  air,  but  the  difference  then  between  a  soil-soil  interface  and  an  air-soil  interface 
should  be  negligible. 
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Berger  and  Doggett  have  shown  that  the  difference  between  a  finite  slab  and 
a  semi-infinite  slab  is  also  very  small.  The  difference  between  the  two  decreases  as 
the  photon  energy  increases  and  as  the  thickness  of  the  slab  increases.  For  example, 
the  dose  build-up  factor  for  a  slab  of  iron  2  mfp  thick  and  for  photon  energy  of  1  Mev 
is  only  10  percent  less  than  for  a  semi-infinite  medium;  for  16  mfp  and  10  Mev  the 
difference  is  only  2  percent. 

The  situation  being  considered  approaches  the  semi-infinite  medium  more  than 
it  does  the  finite  slab  case,  since  the  soil  medium  continues  below  the  shelter  floor 
and  the  minimum  shield  thickness  is  2  mfp.  Consequently,  for  slab  geometry  the 
error  in  using  the  point  isotropic  build-up  factor  is  considered  to  be  less  than  10  percent. 
For  the  other  shelter  shapes  studied,  the  concave  downward  geometry  approaches  the 
semi-infinite  case  even  more  than  slab  geometry  and  is  equal  to  the  semi-infinite  case 
when  the  shelter  radii  equal  zero. 

Since  the  end  result  of  this  report  is  the  geometry  shielding  factors  obtained  by 
the  division  of  two  quantities  containing  these  build-up  factors,  an  inherent  error 
in  them  would  tend  to  cancel  out.  Finally,  Goldstein^  indicates  that  the  accuracy 
of  the  attenuation  coefficient,  fl,  is  even  more  critical  than  the  build-up  factor. 

He  states,  "An  error  of  2%  in  4  means  a  22%  change  in  the  unscattered  dose  at  10  mfp 
and  50%  change  at  20  mfp.  The  uncertainties  in  the  build-up  factors  are  thus  over¬ 
shadowed  by  those  of  jU." 
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BUILD-UP  FAaOR  APPROXIMATIONS 


Goldstein  and  Wilkins  ^  give  a  specific  number  for  fhe  build-up  factors  for 
each  photon  energy,  shield  thickness  (mfp),  and  shield  material.  These  numbers 
can  be  plotted  to  give  curves  of  build-up  versus  energy  or  versus  depth  in  mfp.  In 
order  to  handle  these  functions  as  part  of  the  uncollided  flux  equation  for  gamma 
photons,  the  curves  must  be  represented  by  mathematical  equations.  If  the  build¬ 
up  factor  is  represented  by  an  exponential  equation,  the  manipulation  of  the  basic 
equation  is  simplified  since  it  is  also  exponential  in  nature.  Taylor'^  recognized 
this  v/hen  he  postulated  the  following  equation  for  the  build-up  function: 


Bp  =  A  exp(-a^  fix)  +  (1  -  A)  exp(-a2  fix) 


(7) 


where  A,  a-i,  and  02  are  parameters  dependent  upon  energy  and  material;  and  jix 
is  the  shield  thickness  (mfp).  Table  1  lists  the  values  of  the  parameters  A,  Qi,  and 
02,  which  Taylor  assigned  for  concrete  so  that  this  equation  would  match  the  values 
for  the  point  isotropic  build-up  factors.  Taylor's  equation  will  yield  values  that  are 
within  5  percent  of  the  values  given  by  Goldstein  and  Wilkins. 


Table  I. 


Values  of  A,  and  ^2  for  Concrete  and  a  Point  Isotropic  Source  for  the 
Taylor  Build-up  Function: 

Bp  =  Ae  +  (1  -  A)e 


Photon 

Energy,  Mev 

A 

"°1 

°2 

0.5 

12.5 

0.1110 

1.0 

9.9 

0.0880 

2.0 

6.3 

0.0680 

mmism 

4.0 

3.9 

0.0590 

0.0780 

6.0 

3.1 

0.0585 

0.0830 

8.0 

2.7 

0.0570 

0.0855 

10.0 

2.6 

0.0500 

0.0835 
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Taylor's  equaHon  is  cumbersome,  however,  when  a  numerical  calculation 
must  be  done.  With  over  5,000  separate  calculations  needed  to  obtain  the  results 
of  this  report,  the  use  of  the  Taylor  equation  would  have  increased  this  effort  by  a 
factor  of  two. 

Consequently,  a  simpler  build-up  factor  function  has  been  used: 


B 

r 


(8) 


This  appears  to  be  an  oversimplification,  but,  for  depths  of  2  to  20  mfp,  it  has  proved 
to  furnish  results  within  3  percent  of  results  computed  by  the  Taylor  formula.  The 
values  of  Bq  and  m  were  derived  as  shown  in  Figure  3.  This  figure  indicates  a  build¬ 
up  curve  for  energy  E  versus  depth  in  mfp.  A  straight  line  is  drawn  intersecting  this 
curve  at  2  and  20  mfp.  Bq  then  is  the  intersection  of  this  straight  line  with  the  axis 
and  is  always  greater  than  1.  The  slope  of  the  line  is  m.  Since  Bq  is  a  constant  for 
a  particular  energy  photon,  it  can  be  factored  out  of  the  integral.  The  combination 
of  ju  and  m  results  in  an  "effective"  linear  attenuation  coefficient  jU|.  Values  of  |J.^^ 
versus  energy  are  given  in  Table  II.  Using  this  simplified  build-up  factor  formula, 
all  computations. could  be  done  on  the  basis  of  mean  free  paths  and  all  data  plotted 
as  a  function  of  jif  (mfp).  Consequently,  the  photon  energy  is  not  required  until  an 
actual  problem  is  worked.  The  geometry  shielding  factors  are  valid  then  for  both 
fallout  radiation  or  initial  radiation.  Table  111  is  a  comparison  of  the  geometry  fac¬ 
tor  for  a  sphere  using  the  Taylor  build-up  factor  formula  and  the  simplified  formula 
used  in  this  report.  The  largest  error  is  less  than  3  percent. 

In  view  of  the  fact  that  the  thickness  of  soil  shield  being  considered  will  be 
from  a  jblt  of  2  to  20  mean  free  paths,  the  attenuation  of  the  gamma  photons  in  the 
few  feet  of  air  within  the  shelter  will  not  be  significant  and  will  be  neglected. 

Before  formulating  the  general  problem,  the  following  substitution  will  be 

made: 


H\  =  11-  m 

h 
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Table  II.  Values  of  jl/p,  p.,  p-^/Pt  and  for  Gamma  Phot-on  Energies 
from  1  to  10  Mev  for  Soil  Density  of  1.7  g/cc 


Photon 

p/p  g/cm^ 

iLtft"’ 

2 

p./p  g/cm 

p.  ft  ^ 

1.0 

0.0635 

3.28 

0.0580 

3.00 

2.0 

0.0445 

2.30 

0.0410 

2.12 

4.0 

0.0317 

1.64 

0.0290 

1.50 

6.0 

0.0268 

1.28 

0.0245 

1.27 

8.0 

0.0243 

1.25 

0.0217 

1.12 

10.0 

0.0229 

1.18 

0.0200 

1.03 

Table  III.  Comparison  of  B 


e""  and  B 


E 

pt 

ph 

GF 

GF  Errors 

%  Errors 

1.00 

2 

10 

0.225 

0.220 

-2.2 

1.00 

6 

10 

0.414 

0.425 

+2.1 

1.00 

10- 

10 

0.524 

0.540 

+2.9 

6.00 

2 

10 

0.220 

0.220 

0 

6.00 

6 

10 

0.435 

0.425 

-2.3 

6.00 

10 

10 

0.526 

0.540 

+2.6 

1.00 

2 

15 

0.167 

0.165 

-1.2 

1.00 

2 

20 

0.132 

0.130 

-1.5 
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We  will  now  define  J  as  fhe  integral  portion  of  Equation  2,  The  SfBQ/2  is 
dropped  since  it  is  a  common  factor  for  all  cases.  Thus,  the  general  case  can  be 
represented  by: 


GO- 


I  =  /  '  exp' 

^h  +  t 


9  9  9 

iLt|  t  sec  0  +  h  -|sec  (p  -  (p"^  sin'^  0  cos"^  Q 


m 


-1/2. 

+  sln^  (j)  sln^  0  +  cos^  0)  I 


4L 

R 


(9) 


Figure  4.  Schematic  section  of  shelter  dimensions. 


Using  Figure  4,  It  can  be  shown  that  for  small  Increments  of  d  0  the  following 
identity  holds: 


tan  0  d0 


R 


I 

(10)  i 


11 


I 


Letting  F(0/  0)  =  sin^  (j)  cos^  0  +  sin^  (j)  sin^  0  +  cos^  (j)f  Equation  9 
becomes: 


-1 

exp 


111  tj^  sec^  +  |Ltj  h  sec0  -  u-  h  F  0) 


tan  0  d  (/)  (11) 


Specific  Cases:  Referring  back  to  Figure  2,  the  values  of  t-j  and  ^2 


sec  0 


F(0,  0) 


The  various  specific  cases  can  now  be  generated  by  assigning  specific  values 
to  p  and  q.  The  cases  considered  in  this  report  are: 


(1)  Slab:  p  =  q  =  0 


tn  =  t  sec0 
1  m  ^ 


t2  -  0 


(2)  Sphere:  p  =  q  =  1 


t2  =  h  (sec0  -  1) 
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(3)  Horizontal 
Cylinder: 


p  =  0,  q  -  1 


^2  "  ^ 


9  9  9 

sec0  -  (sin'^^  sin'^0  +  cos'^0) 


-1/2-, 


(4)  Paraboloid: 


p  =  q  7^  1  7^  0 


^2  =  ^ 


2,9 

sec0  -  (q  sin  0  +  cos''"0) 


-1/2, 


Ellipsoid: 


p  /  q  /  1  /  0 


'2  = 


sec0  -  (p^  sin^0  cos^0 


+  q2sin^0sin^0  +  Gos^0) 


-1/2 


(6)  Vertical  Cylinder:  See  Appendix  A. 

SOLUTION  OF  THE  PROBLEM 
Numerical  Integration 

The  solution  of  Equation  11  is  desired  for  the  various  specific  cases  being 
considered.  An  exact  solution  by  integration  is  only  possible  for  the  simple  slab 
and  spherical  cases.  For  the  others,  numerical  integration  must  be  used.  In  order 
to  perform  this  integration  the  source  plane  will  be  divided  into  a  number  of  con¬ 
centric  rings,  centered  over  the  dose  point,  P.  Each  of  these  rings  will  be  determined 


13 


by  the  central  angle  0  to  account  for  the  dependence  of  the  equation  on  0.  The 
plane  will  be  further  divided  into  sectors  to  provide  for  the  dependence  on  the 
angle  0.  This  is  illustrated  by  Figure  5,  Each  small  area  will  then  be  considered 
as  a  separate  dose  source  area^  as  indicated  in  Figure  5,  Area  A^.  This  area  is 
defined  by  the  angles,  0^  -  0n  +  ];  "  ^n  +  1  * 

For  purposes  of  the  numerical  integration,  equal  angle  increments  of  5  degrees 
for  0  were  used  up  to  65  degrees  (contributions  no  longer  significant),  and  15- 
degree  sectors  were  used  for  Q.  The  numerical  integration  was  checked  against  the 
exact  integration  (Table  IV)  for  the  slab  and  sphere  cases  to  test  the  accuracy  of 
using  the  5r:degree  Increments.  The  error  was  less  than  one  percent.  This  is  suffi¬ 
ciently  accurate  for  the  purpose  of  this  study.  Again  the  fact  that  the  slab  case  is 
used  as  a  basis  of  comparison  for  all  other  cases  would  tend  to  cancel  out  any  error 
involved. 


Table  IV.  Comparison  of  Exact  and  Numerical  Integrations  for 
Slab  and  Sphere  Cases 


SLAB 

SPHERE 

Exact 

Numerical 

Exact 

Numerical 

4.89  X  10“^ 
3.78  X  10-3 
3.60  X  10-4 
3.76  X  10-5 
4.16  X  10-5 

4.88  X  10-3 
3.77  X  10-3 

3.58  X  10-4 
3.74  X  10-5 
4.13  X  10-5 

1.04  X  10-3 
1.21  X  10-3 
1.47  X  10-4 
1.76  X  10-5 
2.16  X  lO"*^ 

1.04  X  10-3 
1.21  X  10-3 
1.45  X  10-4 
1.75  X  10“^ 
2.14  X  10-5 
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Figure  5.  Area  division  of  slab  source  for  integration  purposes. 
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Exacf  Integrafion 

The  ini'egraHons  of  the  slab  and  sphere,  based  on  the  same  general  assumptions 
which  have  been  made  for  this  study,  have  been  done  previously  in  NCEL  Report 
TR-025  by  J.  C.  LeDoux,  ^  and  are  merely  presented  below; 


Sf  B 

D  (Slab)  =  —  E,  (H;  tj 


(12) 


D  (Sphere) 


Sf  B 

( 

2 


Mjh 

e  E,  ill.  RJ 


(13) 


where  E^  (Mj  t^)  and  E-j  (/j,]  R^)  are  exponential  integrals  of  the  general  form. 


El  ilix) 


e  ^  dt 
. 


and: 


m 


h  +  t 

m 


These  solutions  were  used  to  check  the  numerical  integrations  for  the  same 
cases  in  order  to  check  the  accuracy  of  the  results  and  the  choice  of  the  5-degree 
increment.  Table  IV  is  a  comparison  of  the  exact  integration  and  the  numerical 
integration.  The  constant  factor  SfBQ/2  has  not  been  included  in  either  case.  The 
difference  in  all  cases  Is  less  than  one  percent.  Based  on  this  comparison,  it  is 
assumed  that  comparable  accuracy  is  maintained  for  all  other  cases. 
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GEOMETRY  FACTOR  CURVES 


The  geometry  factors  (GF)  for  various  shapes,  semi-ellipsoid,  horizontal 
cylinders,  and  vertical  cylinders  are  presented  in  the  form  of  curves. 

Ellipsoids 

Figure  6  is  a  curve  which  fits  the  data  for  the  ellipsoid-type  shelter.  It  has 
a  maximum  error  of  10  percent  and  an  average  error  of  1 .0  percent  from  the  computed 
data  (Appendix  A).  The  curve  is  a  plot  of  GF  vs  t/(t  +  Fh),  where  t  is  the  depth  of 
earth  cover  over  the  crown  of  the  shelter  arch,  h  is  the  height  of  the  shelter,  and 
F  is  an  empirical  factor  which  depends  on  the  ratio  of  the  height  of  the  shelter  to 
the  major  and  minor  axis  radii  of  the  floor  ellipse.  This  function,  t/(t  +  Fh),  was 
used  since  it  produced  the  best  grouping  of  the  ellipsoid  data  on  a  single  curve. 

The  derivation  of  this  function  is  described  in  Appendix  A. 

F  is  obtained  by  entering  the  value  of  (p'^  +  q  ;/2  In  Figure  7. 

where:  p  =  — 

h 

^  =  b 

and;  c  =  the  major  axis  radius 
b  =  the  minor  axis  radius 

The  curve  (Figure  6)  is  also  valid  for  a  sphere  or  a  slab,  since  F  =  1  for  a 
sphere,  and  F  =  0  for  a  slab. 

Horizontal  Cylinders 

Figure  8  is  a  curve  of  GF  vs  t/(t  +  Fh)  for  shelters  of  the  horizontal  cylinder 
type.  This  horizontal  cylinder  does  not  have  to  have  a  uniform  radius.  F  in  this 
case  is: 


since;  q  =  0 

h 

where;  p  =  — 


and:  h  =  height  of  shelter 


c  =  horizontal  distance  from  center  of  shelter  floor  to  outside  arch 
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Geometry  Factor 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

t/(t  +  Fh) 


Figure  6.  Geometry  factors  for  ellipsoid  type  shelters  vs  t/(t  +  Fh), 
where  F  =  [(p2  +  q^)/2  ]  from  Figure  7. 
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Figure  8.  Geometry  factors  for  horizontal  cylinders  vs  t/(t  +  Fh);  where  F  =  p^/2,  and  p  =  h/c. 


Vertical  Cylinders 


Figure  9  is  a  plot  of  the  GF  for  vertical  cylinders  (silos)  vs  r  /h  for  various 
values  of  /ij  f,  where  r^  is  the  radius  of  the  cylinder,  h  is  the  height  of  the  cylinder, 
and  W|  t  is  the  mean  free  paths  of  soil  above  the  roof  of  the  silo.  The  vertical  silo 
case  takes  into  account  radiation  coming  from  the  roof  and  also  from  the  sides. 
Appendix  A  points  out  that  if  the  roof  contribution  alone  is  considered  as  the  contri¬ 
buting  source,  the  maximum  error  would  be  only  5  percent.  Thus,  the  OCDM  manual 
for  aboveground  structures  could  be  used  for  silos  using  only  the  roof  contribution. 

Examples 

Some  examples  of  geometry  factors  for  various  exact  shapes  using  Figures  6, 

8,  and  9  are  given  in  Appendix  B. 


CONCLUSIONS 

The  geometry  factor,  which  has  been  derived  in  this  report,  can  now  be  used 
to  determine  the  additional  nuclear  shielding  provided  by  the  physical  shape  of 
concave  downward  underground  structures  as  compared  to  a  simple  slab  shield.  The 
.  '-“its  of  one  shape  can  be  evaluated  against  that  of  another  and  the  best  one  chosen 
on  the  basis  of  its  shielding,  consistent  with  the  blast  design. 


RECOMMENDATION 

It  is  recommended  that  limited  experiments  be  performed  to  verify  the  results 
of  this  theoretical  study. 
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Appendix  A 
REDUCTION  OF  DATA 


SPHERE,  HORIZONTAL  CYLINDER,  AND  ELLIPSOID  DEVELOPMENTS 

All  of  the  machine  calculations  are  presented  in  Tables  V  and  VI.  This  is 
merely  the  tabulation  of  the  solution  of  Equation  1 1  for  the  various  cases.  This 
data  is  then  used  to  calculate  geometry  factors  (Tables  VII  —  IX)  for  the  various 
cases  considered. 

In  order  to  make  the  results  more  useful,  it  would  be  advantageous  if  the 
geometry  factor  could  be  plotted  as  some  function  of  the  various  parameters  upon 
which  it  depends;  l.e.,  h,  t^,  p,  q,  and  /i|.  The  only  cases  which  can  be  inte¬ 
grated  exactly  ore  the  slab  and  hemisphere.  Consequently,  an  exact  expression 
for  the  geometry  factor  for  a  sphere  can  be  obtained.  Using  Equations  12  and  13, 
we  find  that  the  geometry  factor  for  a  sphere  is: 


GF 


ll^h 

e  El  (mi  Rm) 


El  (M:  f  ) 
1  '^i  m' 


If  we  let:  E]  (x)  =  - F  (x) 

X 


then: 


GF  = 


M|t, 


m 


r  f(mi  R J 


Mi  tm  +  Mi  h  L  H^ll  t^) 


(14) 


*^(^1  '^m)  . 

since:  -  =  1,  the  geometry  factor  for  a  sphere  would  be: 

F(Mi  ^m) 

t 

GF  =  - —  (15) 

(^m  +  h) 


This  is  the  same  geometry  factor  previously  derived  in  Reference  5  for  a  sphere. 
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It  Is  assumed  thai-  the  GF  For  other  shapes  would  have  a  form  similar  to 
Equation  15.  The  actual  attenuation  integral  (Equation  11)  is  very  complex  so 
that  it  is  impossible  to  reduce  the  data  to  a  simple  mathematical  equation  such 
as  Equation  15. 

We  do  know  that  Equation  15  very  closely  satisfies  the  GF  for  a  sphere.  The 
GF  for  a  slab  is  equal  to  one  by  definition.  Equation  15  can  be  modified  to  satisfy 
a  slab  by  the  introduction  of  a  factor  F,  thus: 


GF 


(tm  +  F>') 


(16) 


where:  F  =  0  for  a  slab,  or  1  for  a  sphere 

If  we  examine  Equation  1 1  further  we  notice  that  It  depends  also  on  the  factors 
p^  and  q2.  F  Is  actually  some  perturbation  on  the  height  of  the  shelter  h  and  there¬ 
fore  should  depend  on  the  p^  and  q^  values  in  some  fashion.  Let  us  try: 


p  =  ^ 

2 

This  satisfies  the  slab  and  sphere  conditions  since: 


(17) 


p  =  q  =  0  for  a  slab 
p  =  q  =  1  for  a  sphere 

Using  Equation  17,  a  plot  of  GF  versus  t/(t  +  Fh),  where  t  will  now  be  used 
for  t^,  shows  that  the  computed  data  for  horizontal  cylinders  is  satisfied  by  this 
function  of  F,  since  all  points  plot  very  close  to  a  single  line  (Figure  10).  The 
best  fit  of  these  points  has  been  used  to  produce  Figure  8. 

The  data  for  the  ellipsoid  and  silo  cases  scatters  considerably  from  a  single 
line.  In  ordv3r  to  determine  If  F  could  be  represented  by  some  other  function  besides 
Equation  17,  a  plot  was  made  of  Equation  16  choosing  values  of  F  ranging  from  0.01 
to  10  for  a  juh  of  10,  and  for  various  values  of  4jt.  Figure  1 1  is  this  plot.  Now  if 
the  computed  values  of  GF  from  Table  Vil  are  plotted  on  these  curves,  we  find  that 
they  plot  as  nearly  horizontal  lines.  This  means  that  there  is  some  constant  value 
of  F  which  satisfies  the  various  ellipsoidal  cases. 
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Figure  10.  GF  vs  \/{i  +  F*h)  for  horizonfal  cylinders;  where  F*  =  p^/2  and  p  =  h/< 


9  9 

This  value  of  F  was  then  plotted  against  the  first  guess,  +  q  ;/2,  and 
yielded  a  straight  line  on  log  log  graph  paper,  Figure  7.  Using  the  values  of  F  from 
Figure  7  in  Equation  16  brought  the  values  for  ellipsoidal  shelters  close  to  a  straight 
line.  Figure  12.  The  best  fit  through  these  points  is  the  curve  of  Figure  6.  The  best 
fit  has  a  maximum  error  of  10  percent  and  an  average  error  of  one  percent  from  the 
computed  data. 

Computed  values  of  GF  for  a  jLth  of  15  are  listed  in  Table  Vlll. 


VERTICAL  SILO  DEVELOPMENT 

In  the  case  of  the  vertical  silo,  it  can  be  seen  from  Figure  13  that  t2  =  0 
when  0^0^  and  that  tp  =  h  sec0  -  r^  csc0  when  0  is  greater  than  0^  where 
0(,  =  arc  tan  r^/h . 

Therefore,  the  integration  for  various  values  of  Vq  and  h  was  done  in  two  parts. 
Consider  first  when  <1^  ^  and  t2  =  0.  Equation  11  is  greatly  simplified: 


-00 


1= 


h  +  t 


exp"^  (m  sec0)' 


d  R 


m 


(18) 


where:  R  =  R-j  +  t-j 

R 

from  Figure  13:  sec0 


h  +  t, 


m 


Now  let; 


1  m 


U 


J  h  +  t 


m 


Then  Equation  18  becomes; 


1= 


(•i  +  sec  0 


h  +  t 


m 


d  (II  R) 

exp"'  (il  R)  - J - 

J  ilj  R 


(19) 
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which  is  the  form  of  the  exponential  integral.  It  can  be  easily  shown  then  that: 


11.  (h 


+  t  ) 


-  E. 


4.  (h  +  t  ) 

I  ' 


sec  ( 


(20) 


Substituting  for  the  value  of  4-  in  Equation  20  we  have: 


(11.  tj  -  ill.  t^  sec0) 


(21) 


Now  let  us  consider  when  (j)  >  (j)^.  Equation  11  becomes: 


00 


(h  +  f  )  ^ 


exp 


-flj  (t^  +  h)  sec(p  +  4j  I'o  csc<j) 


d  R 


(22) 


Now  if  we  let  dR/R 


tan  0  d  ^  as  before,  then  Equation  22  becomes: 


“^i  (’’m  ■*'  ^^i  ‘'o 


tan  ^  d  (/) 


(23) 


Equation  23  was  numerically  integrated  by  machine,  using  increments  for  0  of  5  degrees 
until  the  contribution  was  negligible.  The  sum  of  the  results  from  Equation  23  and 
Equation  21  for  each  case  was  divided  by  the  slab  case  E-]  (fij  t^)  in  order  to  give 
the  geometry  factors  as  defined  above. 
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Figure  13.  Schematic  of  vertical  silo  configuration. 


Table  IX  lists  the  geometry  factors  obtained  for  various  values  of  fQ,  h^  and 
I'ffl*  Figure  9  graphically  displays  the  resultS;^  plotting  geometry  factor  versus 
the  ratio  Tq/U  for  various  thicknesses  of  material  in  mfp.  The  last  column  of  Table  IX 
is  the  geometry  factor  when  only  Equation  21  is  used.  It  can  be  seen  that  the  largest 
variation  is  only  5  percent,  occurring  at  one  of  the  less  important  cases. 

The  value  of  |l'^  t^  is  calculated  for  the  particular  energy  involved.  Table  II 
lists  the  values  of  ll\/p,  and  flj  for  soil  of  density  1.7  gm/cc  versus  photon  energy. 
Knowing  the  density  of  the  soil  to  be  used  and  multiplying  /ij/p  by  this  density, 

Pj  t^  can  be  calculated  in  any  appropriate  units  for  use  as  input  to  Figure  9. 
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XVlIl  0.0488  0.0819  0.1119  0.1394  0.1650 

XIX  0.0332  0.0563  0.0775  0.0974  0.1161 

XX  0.0224  0.0384  0.0533  0.0676  0.0812 
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Appendix  B 

EXAMPLES  OF  CO/liJ°UrED  GEOMETRY  FACTORS  FOR  VARIOUS  SHELTERS 


1.  Horizonlal  Cylinder 


f 

r  +  Fh 


5 

- ij -  =  0.455 

5  +  i(12) 


2.  Sphere 


Thus,  from  Figure  6,  GF  =  0.348 
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3.  Ellipsoid 


12 

—  =  0,8 

15 


U 

25 


=  0.48 


p2  +  q2  0.64  +  0.23 


From  Figure  7,  F  =  0.475 


t  ^  5  _ 

f  +  Fh  5  +  (0.475)  12 
Thus,  from  Figure  6,  GF  = 


4.  Caffle  Pass 
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—  =  2.677 
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—  =  3.56 

2 

From  Figure  7,  F  =  3 


f  +  Fh  5  +  3(8) 

From  Figure  6,  GF  =  0.22i 


=  1.27  (5)  =  6.35 
From  Figure  9,  GF  =  0.565 


1 

I 

I 

i 


43 


DISTRIBUTION  LIST 


No.  of 

SNDL 

Copies 

Code 

10 

Chief,  Bureau  of  Yards  and  Docks 

BuDocks  Standard  Distribut'-m 

1 

23A 

Naval  Forces  Commanders  (Taivan  Only) 

2 

39B 

Construction  Battalions 

9 

39D 

Mobile  Construction  Battalions 

3 

39  E 

Amphibious  Construction  Battalions 

2 

39  F 

Construction  Battalion  Base  Units 

1 

A2A 

Chief  of  Naval  Research  -  Only 

2 

A3 

Chief  of  Naval  Operations  (Op-07,  Op-04) 

5 

A5 

Bureaus 

3 

B3 

Colleges 

2 

F4 

Laboratory  ONR  (Washincton,  D.  C.  Only) 

1 

E16 

Training  Device  Center 

8 

F9 

Station  -  CNO  (Boston;  Key  West;  New  Oiteans;  Son  Juan;  Long  Beach;  Son  Diego; 
Treasure  Island;  and  Rodman,  C.  Z.  Only) 

5 

F17 

Communication  Station  (Stin  Jiion;  Son  Francisco;  Pearl  Horbor;  Adah,  Alaska;  and 
Guam  only) 

1 

F21 

Administration  Command  and  Unit  CNO  (Saipon  only) 

2 

F40 

Communication  Facility  (Pt  L/autey  and  Jopon  only) 

1 

F41 

Security  Stotion 

2 

F42 

Radio  Station  (Oso  and  Choltanhom  only) 

1 

F48 

Security  Group  Activities  (Win  er  Horbor  only) 

8 

H3 

Hospital  (Chelsea;  St.  Albon'i;  Portsmouth,  Vo;  Beaufort;  Groat  Lakes;  Son  Diego; 
Oakland;  and  Camp  Pendleton  only) 

1 

H6 

Medical  Center 

2 

J1 

Administration  Command  and  Unit-BuPers  (Groot  Lokes  and  Son  Diego  only) 

1 

J3 

U.  S.  Fleet  Anti-Air  Worfore  Tr  ^'ining  Center  (Virginia  Beach  only) 

2 

J4 

Amphibious  Bases 

1 

J19 

Receiving  Station  (Brooklyn  only) 

1 

J34 

Station  -  BuPers  (Washington,  D.  C.  only) 

1 

J37 

Training  Center  (Bainbridge  only) 

1 

J46 

Personnel  Center 

1 

J48 

Construction  Training  Unit 

1 

J60 

School  Academy 

/ 

i 


44 


Distribution  List  (Cont*d) 


No.  of 

SNDL 

copies 

Code 

1 

J65 

School  C£C  Officers 

1 

j84 

School  Postgraduate 

1 

J9C) 

School  Supply  Corps 

1 

J95 

School  V/ar  CoMege 

1 

j99 

Communi cation  Training  Center 

11 

LI 

Shipyards 

4 

L7 

Laboratory  •  3uSh'ps  (New  London;  Panama  Lity;  Carderock;  and  Annapolis  only) 

5 

L26 

Naval  Fociiities  -  BuShlps  (Antigua;  Turks  Island;  Barbados;  San  Salvador;  ond 
Eleuthera  only) 

1 

L30 

Submarine  Bo  -.e  (Grolon,  Conn,  only) 

2 

L32 

Naval  Support  Activities  (London  &  Noples  only) 

2 

L42 

Fleet  Activities  -  BuShlps 

4 

M27 

Supply  Center 

7 

M28 

Supply  Depot  (Except  Guantanamo  Boy;  Subic  Boy;  and  Yokosuka) 

2 

M61 

Aviotion  Supply  Office 

3 

N1 

BuDocks  Director,  Overseas  Division 

42 

N2 

Public  Works  Offices 

7 

NS 

Construction  Bottolion  Center 

5 

N6 

Construction  Officer-in-Chorge 

1 

N7 

Construction  Rr  sident-Officer-in-Chorge 

12 

N9 

Public  Works  Center 

1 

N14 

Housing  Activi* 

2 

R9 

Recru’t 

•-rr-/  ;uns  ^^loany  and  uurstow  only;  ▼  ▼  ▼  V 

1 

R20 

Marine  Corps  Schools,  Quantico 

3 

R64 

Marine  Corps  B'^se 

1 

R66 

Marine  Corps  Camp  Detachment  (Tengan  only) 

7 

WlAl 

Air  Station 

35 

W1A2 

Air  Station 

9 

WIB 

Air  Station  Auxiliary 

5 

WIC 

Air  Facility  (Pi  oenix;  Monterey;  Oppomo;  Noho;  and  Naples  only) 

3 

WIE 

Marine  Corps  Air  Stotion  (Except  Quantico) 

1 

V/IF 

Marine  Corps  Auxiliary  Air  Station 

8 

WIH 

Station  -  BuWeas  (Except  Rota) 

Distribution  List  (Cont’d) 


Mo.  of 
oop  ie  s 

1 

1 

I 

1 

1 

1 

I 

1 

10 

1 

3 

2 

2 

2 

10 

2 

2 

2 

2 

1 


Chief  of  Staff,  U.  S.  Army,  Chief  of  Research  and  Development,  D<  part.Tienf  of  the  Army, 

Washington  2.5,  D.  C. 

Office  of  the  Chief  of  Engineers,  Asst.  Chief  of  Engineering  for  Ci-  ,1  Vforks,  Department  of  the 
■Army,  Washington  25,  D.  C. 

Chief  of  Engineers,  Department  of  the  Army,  Attn;  Engineering  R  &  D  Division,  Washington  25,  D.  C. 

Commanding  Officer,  Engineering  R  &  D  Laboratories,  Attn:  Techn'cal  Intelligence  Branch, 

Fort  Selvoir,  Virginia 

Commanding  General,  Wright  Air  Development  Center,  Air  Research  ..  id  Development  Command. 

Wright- Patterson  Air  Force  Base,  Ohio 

Deputy  Chief  of  Staff,  Development,  Director  of  Research  and  Developm-. ,nt.  Department  of  the 
Air  Force,  Washington 

President,  Marine  Corps  Equipment  Board,  Marine  Corps  Schools,  Quoni.:o,  Virginia 

Director,  National  Bureau  of  Standards,  Department  of  Commerce,  Corun  :ticut  Avenue, 

Washington,  D.  C. 

Armed  Services  Technical  Information  Agency,  Arlington  Hall  Station,  Arlington  12,  Virginia 

Deputy  Chief  of  Staff,  Research  and  Development  Headquarters,  U.  S  tarine  Corps 

Headquarters,  USAF,  Directorate  of  Civil  Engineering,  Attn;  AFOCE  !  S,  V/oshington  25,  D.  C. 

Commander,  Hsadquarters,  Air  Peseorch  and  Development  Command,  A  nlrevvs  Air  Force  Base, 
Washington  25,  D.  C. 

Office  of  the  Director,  U.  S.  Coast  and  Geodetic  Survey,  V/ashington  25,  D.  C. 

Library  of  Congress,  V/ashington  25,  D.  C. 

Director,  Office  of  Technical  Services,  Department  of  Commerce,  V/ashir  gton  25,  D.  C. 

NCEL  Standard  Distribution 

Director  of  Defense  Research  and  Engineering,  Deportment  of  Defense,  V-ashington  25,  D.  C. 

Director,  Division  of  Plans  ond  Policies,  HoodqMorters,  U  S.  Morinc  Corps,  Woshington  25,  D,  C. 

.  -t. . . . 

Commanding  Officer,  U.  S.  Nava!  Construction  Battalion  Center,  Attn!  Tcchnicol  Division, 

Code  141,  Port  Hueneme,  California 

Commanding  Officer,  U.  S.  Navol  Construction  Battalion  Center,  Attn:  Materiel  Department, 

Code  142,  Port  Hueneme,  Californio  * 

Commanding  Officer  (Patent  Dept.),  Office  of  Naval  Research  Branch  Office,  1030  E.  Green  Street 
Pasadena,  California  * 


46 


Distribution  List  (Cont’d) 


NCEL  Supplemental  Distribution 

Commandant,  Industrial  College  of  the  Armed  Forces,  Washington,  D.  C. 

Commandant,  U.  S.  Armed  Forces  Staff  College,  U.  S.  Naval  Base,  Norfolk,  Va. 

Chief,  Bureau  of  Ships,  Attn:  Chief  of  Research  and  Development  Division,  Navy  Deportment 
Washington,  D.  C. 

Officer  in  Charge,  U.  S.  No'/y  Unit,  Rensselaer  Polytechnic  Institute,  Troy,  N.  Y. 

Chief  Bureau  of  Naval  Weapons,  Attn:  Research  Division,  Navy  Department,  Washington,  D  C 
Commanaer,  Pacific  Missile  Range,  Attn:  Technical  Director,  Point  Mugu,  Calif. 

Commanding  Officer,  U.  S  Naval  Unit,  U.  S.  Army  Chemicol  Corps  School,  Fort  McClellan  Ala 
Deputy  Chief  of  Staff,  Research  -S.  Development  Headquorters,  U.  S.  Marine  Corps,  Washington  0  C 
Deputy  CCMLO  for  Scientific  Activities,  Washington,  D  .C. 

Chief  of  Ordnance,  U.  S.  Army,  Attn:  Research  &  Development  Loborotory,  Washiriuij',.  ^ 

U.  S.  Army,  Attn:  Director  of  Research  and  Development  Group,  Washington.  C. 

Directorate  of  Medical  Research,  Chemical  Warfare  Board,  Army  Chn,ilcal  Cente.',  Md. 

U.  S.  Aimy  Corps  of  Engineers,  Office  of  the  District  unui-  ..er,  S).  poul  District,  1217  U  S  P  0 
and  Customs  House,  St.  Paul,  Minn. 

Snow,  Ice,  and  Permafrost  Research  Fv,tnL!i"hmpnt,  Corps  of  Engineers,  U.  S.  Army,  l21S  Washington 
Avenue,  Willmette,  III. 

Taft  Sanitary  EnginecMrtn  Csnrer  U5!^*HS,  4676  Columbia  Parkway,  Cincinnati,  Ohio 

Commander,  Air  Rese-’iLn  &  Development  Commond,  Attn:  Library,  Andrews  Air  Force  Bose 
Washingioii,  D. 

Directorate  of  Research,  Air  Force  Special  Weapons  Center,  Kirtlond  Air  Force  Bose,  N.  M 
Commanding  Officer,  Biologicol  Warfare  Laboratories,  Fort  Detrick,  Frederick,  Md. 

Ssndio  CorpoiOtion,  Attn-  Classified  Document  Division.  Box  SftOn^  ihuquc  rqin  ,  h.  jyj 
--c-^i^vafTdvct^n^  Oniversny'C^oiii^^  «ic,n..ond  |-ieT» Station,  Blakeley,  Calif. 

Library,  Engineering  Department,  Stanford  University,  Stonford,  Colif. 

Library,  Harvard  University,  Graduate  School  of  Engineering,  Combridge,  Mass. 

Director,  Engineering  Research  Institute,  University  of  Michigan,  Ann  Arbor,  Mich. 

Library,  Engineering  Department,  University  of  Californio,  405  Hiigard  Avenue,  Los  Apgeles 
Library,  University  of  Southern  Californio,  University  Pork  Los  Angeles 
Director,  Marine  Laboratory,  University  of  Miami,  Coral  Gables,  Fla. 

Director,  Soil  Physics  Laboratory,  Department  of  Engineering,  Attn:  Library,  Princeton  University, 
Princeton,  N.  J. 

Chief,  Bureau  of  Yards  and  Docks,  Code  0-230,  Washington,  0.  C. 

Mr.  W.  R.  Perret,  5112,  Sandia  Corporation,  Sandio  Bose,  Albuquerque,  N.  M. 

Dr.  N.  M.  Newmark,  Civil  Engineering  Hall,  University  of  Illinois,  Urbana,  111. 

Mr.  Fred  Sauer,  Physics  Department,  Stanford  Research  Institute,  Menlo  Park,  Calif. 

Dr.  Harold  Brode,  RAND  Corporation,  1700  Main  Street,  Santa  Monica,  Calif. 


47 


Distribution  List  (Cont’d) 


Mr.  Kenneth  Kap'on,  Broodvievr  i^eseorch  Corporation,  1811  Trousdale  Drive,  Burlingame,  Calif. 

Prof.  J.  Neils  Thompson,  Civil  Engineering  Department,  University  of  Texas,  Austin,  Tex. 

Mr.  G.  L.  Arbuthnot,  Waterways  Experiment  Station,  Post  Office  Box  631,  Vicksburg,  Miss. 

Mr.  William  J.  Taylor,  Terminal  Ballistics  Laboratory,  Aberdeen  Proving  Ground,  Aberdeen  Proving 
Ground,  Md. 

Dr.  T.  H.  Schiffmon,  Armour  Research  Foundation  of  Illinois  Institute  of  Technology,  Technology 
Center,  Chicago,  III. 

Chiefof  Engineers,  Deportment  of  the  Army,  Washington,  D.  C. 

Mr.  Eric  Wang,  Air  Force  Special  Weopons  Center,  Kirtlond  AFB,  Albuquerque,  N.  M. 

Commander,  Air  Force  Ballistic  Missile  Divsion,  Air  Research  and  Development  Command, 

Attn:  Or.  George  Young,  Post  Office  Box  262,  Inglewood,  Calif. 

Dr.  James  0.  Buchanan,  Office  of  Civil  and  Defense  Mobilization,  Bottle  Creek,  Mich. 

Mr.  L.  Neal  FitzSimons,  Office  of  Civil  and  Defense  Mobilization,  Winder  Building,  Washington,  D.  C 

Mr.  E.  E.  Sholowitz,  Protective  Construction,  GSA  Building,  19th  and  F  Street,  N.  W,, 

Washington,  D.  C. 

Lt.  Col.  Russell  J.  Hutchinson,  052921,  Office  of  Area  Engineer,  Saudi  Arabia,  U.  S.  A,  Engineer 
District,  Trans-East,  APO  616,  New  York 

Capt.  A.  B.  Chilton,  CEC,  USN,  U.  S.  Naval  Civil  Engineering  Laboratory,  Port  Hueneme,  Calif. 

Cdr.  C.  A.  Grubb,  CEC,  USN,  Public  Works  Center,  Navy  No.  128,  FPO,  San  Francisco 

W.  M.  McLollono,  CEC,  USN,  Public  Works  Office,  U.  S,  Naval  Base,  Charleston,  S,  C. 

Copt  L.  N,  SoLMders,  Jr.,  CEC,  USN,  Bureau  of  Yords  onci  Docks,  Code  D-400,  Washington,  D.  Ct 
Cdr.  J,  F,  CIOtKC,  CEC,  ubN,  Area  Public  Works  Office,  Chesapeake,  U.  S.  Naval  WenDons  Plan* 

'  * .  I  r*  <*• 

~  .,,,1  ../istlii..,  isavy  i-io.  iZo, 

FPO  San  Fro, 1C  ■,  a 

Mr.  G.  H.  AlbriqK^.  The  Pcnnsylvonio  State  University,  College  of  Engineering  and  Architecture, 
University  Park,  i’ci'i, 

Capt.  J,  H.  Barker,  Jr,  CEC.  USN,  U.  S.  Naval  Missile  Center,  Point  Mugu,  Calif. 

Cdr.  W.  J.  Christen  .cn  CEC,  USN,  Bureau  of  Yards  and  Docks,  Code  D-440,  Washington,  D.  C. 

Capt.  J.  H.  Lofland,  Jr.,  CEC,  USN,  Pearl  Harbor  Naval  Shipyard,  Novy  No.  128,  FPO, 

San  Francisco 

Capt.  W.  A.  McManus,  CEC,  USN,  U.  S.  Naval  Air  Station,  Norfolk,  Vo. 

Lcdr.  J.  D.  Androv/s,  CEC,  USN,  Defense  Atomic  Support  Agency,  Washington,  D.  C. 

Cdr.  0.  P.  Cunning,  CEC,  USN,  District  Public  Works  Office,  4th  Naval  District,  Nava!  Base, 
Philadelphia,  Penn. 

Lcdr.  V/.  J.  Fraucy,  CEC,  USN,  Bi''e<  o  of  Yards  and  Docks,  Code  D-500,  Washington,  D.  C. 

Lcdr.  F.  W.  Galbraith,  CEC,  USN,  U  S.  Naval  Air  Stotion,  North  Island,  Son  Diego 


48 


